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, This mini-review discusses the recent contribution of theoretical and computational physics as well as experi- 

W ■ mental efforts to the understanding of the behavior of colloidal particles in confined geometries and at liquid 

^ ' crystalline interfaces. Theoretical approaches used to study trapping, long- and short-range interactions, and 

assembly of solid particles and liquid inclusions are outlined. As an example, an interaction of a spherical 
colloidal particle with a nematic-isotropic interface and a pair interaction potential between two colloids at 
this interface are obtained by minimizing the Landau-de Gennes free energy functional using the finite-element 
method with adaptive meshes. 
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1. Introduction 

>, 

Liquid crystals are known for their anisotropic mechanical and optical properties which originate 
from the long-range orientational molecular ordering. If a liquid crystal is used as a host liquid 
in a colloidal suspension, this ordering gives rise to additional long-range interactions between 
the colloidal particles. The type of the interaction is controlled by the presence and symmetry of 
topological defects of the director field. Particle clustering, formation of superstructures, and even 
new phases are immediate consequences of these anisotropic interactions. 

In what follows, we present a brief review of theoretical and simulation methods which are 
used to study liquid crystalline systems. Each description of phenomena occurring on a particular 
length-scale has its own order parameter, such as the director field, tensorial order parameter, or 
particle density. For example, long-range forces between colloidal particles can be calculated by 
direct integration over the director field around them. For smaller particle-particle separations, 
where nonlinear effects become important and topological defects change their positions, a direct 
minimization of the Landau-de Gennes free energy functional with respect to the tensorial order 
parameter can be used. For even smaller separations, presmectic ordering and density variations 
cannot be ignored and density functional approaches or computer simulations can be employed. 

With these techniques at hand, we discuss different physical situations: (i) a single particle 
suspended in a nematic host; (ii) an effective pair interaction of two particles suspended in a 
bulk nematic; (iii) an interaction of a particle with a nematic-isotropic (NI) interface and (iv) 
an interaction of two colloidal particles captured by this interface. We also give an overview of 
colloidal self- and directed-assembly in two-dimensional systems, such as free-standing smectic films 
and quasi two-dimensional nematic colloids. Finally, to give an example, we calculate the expelling 
force on a spherical particle in the vicinity of the NI interface and the effective interaction potential 
between two spherical particles captured by this interface. In both cases we assume the film-like 
shape of the domain occupied by the nematic phase and study how the film thickness influences 
effective interactions. 
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2. Theoretical background 

2.1. The order tensor and the nematic director 

To form even the simplest liquid crystalline state one needs anisotropic molecules. To a good 
approximation, these highly anisotropic molecules can be thought of as rigid rods or ellipsoids of 
revolution. The orientation of every single molecule is then fully described by a unit vector u along 
its anisotropy axis. 

Given the orientations of all the molecules one can define an orientational order tensor 



S a p (r) = — ]T 



(1) 



where the sum is over the all N molecules in a small but macroscopic volume located at the point 
r, itW is a unit vector along the axis of the z-th molecule. The order tensor S is a symmetric 
traceless tensor and in general case has five independent components. Since the anisotropy of a 
macroscopic state originates from the molecular anisotropy on a microscopic level, the orientational 
order tensor Eq. ((T|) determines the anisotropy of macroscopic properties of mesophase. 

The director n is an eigenvector of the order tensor S a p corresponding to the largest eigenvalue 
S. It gives an average orientation of the molecules in a mesoscopic volume located at the point 
r. The order parameter S has clear physical meaning: in the isotropic liquid S = 0, a perfectly 
aligned system has order parameter 5 = 1, and S = —1/2 if the molecules are randomly aligned 
in a plane. 



2.2. Frank-Oseen elastic free energy 

In the case where the scalar order parameter S is constant, an expansion in terms of the director 
n is normally used to calculate elastic free energy density of a nematic [ [1] 

gt = i [k x (V • nf + K 2 (n • [V x n]f + K 3 [n x [V x n]] 2 } , (2) 

where K\, K 2 , and K 3 are Frank elastic constants representing splay, twist, and bend director 
distortions, respectively. 

Equations for the director can be obtained by minimization of the free energy functional J v gtdr 
with an additional constraint n 2 — 1 and have the following form 

- KxV (V • n) + K 2 {A[V x n] + [V x (An)}} + K 3 {[[V x n] x B] + [V x [B x n]]} + /m = 0, (3) 

where A = n ■ [V X n], B = [n x [V X n]], /x is an unknown Lagrange multiplier which must be 
found from the condition n 2 = 1. 



2.3. Landau-de Gennes free energy 

Following the general approach to phase transitions given by Landau theory [[2] one can assume 
that the Gibbs free energy density g (P, T, S a p) is an analytic function of the order parameter S a p. 
Expanding g in a power series and taking into account invariance of g under rotations one can 
obtain [[3] 

ffbulk = aTrS 2 - bTvS 3 + c [TrS 2 ] 2 , (4) 

which is the most general expansion correct up to the fourth order in S a p- The coefficients a, 6, 
and c are in the general case functions of pressure P and temperature T. In practice, a is assumed 
to depend linearly on the temperature, b and c are considered temperature independent. 

For this free energy, the uniaxial nematic state is stable when b 2 > 2Aca with the degree of 
orientational ordering 

b / r 64^\ 
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The order tensor S can depend on the spatial coordinates, which means that either the director 
n or the scalar order parameter S (or both) vary from point to point. This variation can be due 
to external forces imposed on the system, thermal fluctuations, or boundary conditions. Assuming 
these deformations to vary slowly in space relative to the molecular distance scale, it is possible 
to describe the response of the liquid crystal using continuum elastic theory. Then the Gibbs free 
energy density, relative to the free energy density of the state with the uniform orientational order 
can be written as 

dSij dSij dSij dSik , T dSij dSik , . 

gd = L 1 —^-—^-+L 2 —^-- — + L 3 — , (6) 

OXk OXk OXj OXk OXk OXj 

where L\, L 2 , and L 3 are the elastic constants and the summation convention is assumed. 



2.4. Density functional theory 

Onsager theory [ [4] is a form of a molecular field theory, with a distinctive feature that the 
configurational average energy is replaced by the orientation-dependent configurational entropy. It 
is a forerunner of the modern density- functional theory [ [5] . 

Consider an ensemble of elongated particles interacting pairwise through some potential V(l, 2) 
which depends on both positions and orientations of the molecules 1 and 2, where 1 = (t*i, fii) and 
2 = (t"2, $12). Onsager has shown how the Mayer cluster theory may be used to give an expansion 
for the equation of state of this system [3]. Onsager's expression for the Helmholtz free energy 
per particle is expressed in terms of the single-particle density, p(l), which gives the number of 
molecules per unit solid angle and per unit volume 

PF\p) = J p(l) {lnp(l)A 3 - 1 - An + (3U(1)} d(l) - l -J f(l, 2)p(l)p(2)d(l)d(2). (7) 

Here /3 = with fee being the Boltzmann constant and T the absolute temperature, A is the 

de Broglie wavelength, fi is the chemical potential, U(l) is the external potential energy, 

/(l,2)=exp[-V(l,2)/fc B T]-l (8) 

is the Mayer /-function, and V(l, 2) is a pair-wise interaction potential. 

The equilibrium single-particle density that minimizes the free energy functional in Eq. ([7]) is 
the solution of the following Euler-Lagrange equation 

\np(l)A 3 -f3p + (3U(l)- J /(l,2)p(2)d(2) = 0, (9) 

which can be obtained from the variation of the functional ([7} . 

In principle, Eq. ([9]) can be solved numerically for any form of the second virial coefficient which 
can be expanded in Legendre polynomials. However, sometimes it is more convenient to minimize 
the functional in Eq. (0 directly instead of solving the integral equation (j9]). The solutions are 
qualitatively the same as the solutions of the Maier-Saupe equation [|3]. The difference is that the 
Maier-Saupe theory is commonly applied to liquids, which are only slightly compressible, whereas 
the Onsager expansion is applied to dilute suspensions of particles for which the change in free 
energy with density is relatively small. The Onsager theory predicts large changes in density at 
the transition, as is observed for such systems. Note that the Onsager theory does not describe the 
bulk equation of state perfectly. The reason for this is the truncation of the virial expansion of the 
free energy at the leading, pairwise term. However, systematic improvements are possible which 
lead to the better agreement with the bulk equation of state [ [SJ [7J |H1 IS] ■ 



2.5. Monte Carlo and molecular dynamics simulations 

Computer simulations is a well established method to study liquid crystalline mesophases. It 
has been used to simulate various bulk properties of liquid crystals, such as bulk elastic constants [ 
□33 BU [HI Q21 El EH El, viscosities [EHH, helical twisting power [ [H [20], parameters of the 
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(a) (b) (c) (d) 

Figure 1. (Color online) Sketches of (a) the satellite defect (dipole) and (b) the Saturn-ring 
defect (quadrupole) in the vicinity of a sphere with normal surface anchoring placed in a uniform 
aligned nematic. (c) and (d) show the corresponding scalar order parameter profiles and director 
maps obtained from the molecular dynamics simulations. Adapted from Ref. [ 31. 

isotropic-nematic interface [ 21 . Moreover, recent works have also proved that computer simulations 
can be effectively applied to study confined geometries [ [22] [23] . 

Computer simulations can also be used to study topological defects in liquid crystals. They 
provide a direct access to molecular positions and orientations and hence can directly validate 
aforementioned macroscopic and mesoscopic descriptions. In many cases, computer simulation 
provides information not extractable from phenomenological theories, such as the exact structure 
of the disclination line core [[24j 25 or about the detailed liquid crystal structure around external 
inclusions [ [26j [27] , defects structure and resulting effective interactions in nematic colloids [ 28, 
1291 130] , confirming the corresponding predictions of the continuous Landau-de Gennes theory [ 
EH [29]. 

3. An isolated colloidal particle 

In order to understand an effective interaction between particles dispersed in a nematic liquid 
crystal, it is essential to know the liquid crystal ordering near one of the particles. An isolated col- 
loidal particle can provide either homeotropic (perpendicular to the surface) or tangential (parallel 
to it) boundary conditions for the nematic director. A particle with sufficiently strong homeotropic 
anchoring carries a topological charge of strength +1. If the director field is uniform far away from 
the particle, i. e. the total charge of the whole system is zero, an additional defect in the medium 
must be created in order to compensate the topological charge of the particle. Two defect types 
can arise in this case. The first one is a hyperbolic hedgehog, a point defect characterized by the 
topological charge of —1, called satellite defect. This particle-defect pair is usually modeled as a 
topological dipole because the asymptotic behavior of the director field has dipolar symmetry. The 
other possibility is a quadrupolar, or a Saturn-ring defect, that is a —1/2 strength disclination ring 
that encircles the particle. Both are sketched in Fig. [T] 

Theoretical and numerical work based on the Frank- Oseen elastic free energy [ [32] [33] [34] [35] 
as well as experimental observations [ [36] suggest that the dipolar configuration is stable for the 
micron-sized particles usually considered experimentally. The Saturn-ring configuration should 
appear if the particle size is reduced and, when present, is always predicted to be most stable in 
the equatorial plane normal to the far field director. If the strength of the surface anchoring is weak 
enough, a surface- ring director configuration is possible [ [37] [34] . It has also been found that a 
twisted dipolar configuration is possible for high enough ratios of splay to twist elastic constants [ 
[34]. 

In order to describe the defect structure, the elastic theory based on either Frank-Oseen, Eq. ([2]), 
or Landau-de Gennes, Eqs. (j4j), ([6]), expressions for the free energy can be used. Note that the Frank- 
Oseen elastic theory has several restrictions: the defect core region cannot be treated and the liquid 
crystalline phase is considered as uniaxial. It has been demonstrated [ 38, 39, 40, 24, 26 that the 
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defect core of the simplest disclination line of strength +1 has a rather complex structure: the order 
parameter varies near the core region, the uniaxial phase becomes biaxial, and even the density of 
the mesophase can be affected by the presence of a wall or a defect [ [3TJ [22] . Hence, for a better 
description, it is necessary to use the tensorial order parameter formalism [ [41] |42j [43] |44j [45] . This 
description has no anomalous behavior at the singularity of the director field and can resolve order 
parameter variations as well as biaxiality of the defect core. In addition, density variations can be 
taken into account using Onsager theory [ [35] . 

4. Effective pair interactions 

Several types of interactions exist between the particles in a colloidal suspension: van der Waals, 
screened Coulombic, steric, fluctuation-induced Casimir forces, depletion, etc. In a nematic solvent, 
elastic deformations of the director around particles lead to an additional, long-range, interaction. 
It can be of a dipolar or quadrupolar type, depending on the symmetry of the director configuration 
around particles. 

The first experimental observations of these interactions were made in inverted nematic emul- 
sions, that is water droplets dispersed in a nematic solvent [US] - The colloidal droplets formed linear 
chains, which were breaking upon transition to the isotropic phase. It was concluded that the liquid 
crystal host induces dipolar interaction which leads to droplet chaining, as well as a shorter range 
repulsive interaction [ 147). Later, both dipolar and quadrupolar colloidal forces were measured 
directly using ferrofluid droplets in a magnetic field [HE], magneto-optic tweezers [ [49] [50] [51] , iron 
particles in a magnetic field [ , and dual-beam optical tweezers [ [53] [M] [S3] . 

We will first discuss the asymptotic behavior (large inter-particle separations) of forces mediated 
by a nematic host. Possible corrections at small separations due to reorientation of defects and 
higher order terms in the expansion of the elastic free energy are presented at the end of this 
section. 



4.1. Long-range asymptotic behavior 

The large distance behavior of nematic-mediated interactions can be obtained from an electro- 
static analogy by assuming that all constants in the expression for the elastic free energy, Eq. , 
are equal (one-constant approximation, K\ = K2 = Kg = K) and by performing a multipole 
expansion of the director field [ 33 5 56! ■ Indeed, far away from a colloidal particle, the director can 
be written in a linearized form, n(r) = (n x ,n y ,l — 0(n^.,nt)), and the Frank-Oseen elastic free 
energy, Eq. ([2]), reads as 

d 3 r((Vn x ) 2 + (Vn,) 2 + 0«,^)). (10) 

Hence, the transverse director components, n x and n y , fulfill the Laplace equation. For a single 
particle, these can be expanded into multipoles. Using the superposition approximation, one can 
derive the following expressions for the effective pair potential 

1 - 3 cos 2 e 

•"aipolc— dipolc OC g , 

9 -90 cos 2 6» + 105 cos 4 9 

^quadrupolc — quadrupolc ^ ^5 7 (1 1) 

COS 9 2 
Vdipolc-quadrupolc °C ^4 (15 COS 9 — 9), 

where d is the distance between the particles and 9 is the angle between the far field orientation of 
the nematic director (the z axis in our case) and the vector connecting the centers of the colloidal 
particles. These expressions can be generalized for particles of arbitrary shapes [ [57] [551 [52] or for 
weak anchoring at a particle surface [ [60] . 
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The asymptotic behavior of the effective interaction potentials, given by Eqs. (fTTj) . has been 
directly measured for the dipole-dipole configuration corresponding to 9 = 0, using the free-release 
and the dual-beam optical tweezing methods. For parallel dipoles, this force was found to be 
attractive and decaying with the particle-particle separation as d~ 4 [ [54], [53 [61} [53] . However, for 
the antiparallel orientation of dipoles, the pre-asymptotic decay was observed and the repulsive 
force decayed as d~ 36 [ 55] . Numerical minimization of the Landau-de Gennes free energy functional 
predicted an even slower decay of d~ 3 [ I62U63) . In this work, the force was evaluated by numerically 
differentiating the Landau-de Gennes free energy. The direct integration of the corresponding stress 
tensor resulted in a decay of d~ 363 [[64]- Although the distance between the particles was as large 
as five particle diameters, the anti-parallel alignment was still in a pre-asymptotic regime. The 
dipole-quadrupole attraction with a decay of d~ 5 , corresponding to 9 — it in Eqs. (jlip . was also 
confirmed experimentally [ 65 a . The quadrupolar tail of d~ e was also found experimentally for 
colloidal particles with the tangential boundary conditions [[2§JGI2]- Theoretically, the quadrupole- 
quadrupole asymptotic behavior was studied in detail using molecular theory, which accounts for 
the anisotropy of the molecular interactions [ 67 . 



4.2. Short-range behavior 

The electrostatic analogy of the linearized Frank-Oseen free energy in the one-elastic-constant 
approximation, Eq. (1101) , gives a good physical insight to understanding of structural properties of 
particle assemblies in dipolar [[61], quadrupolar [[68], or mixed [[65] [69] situations. For example, 
the formation of long chains in the first, or rhomboidal structures in the last two cases can be 
inferred from the 0-angle dependence of the potentials in Eqs. (| 1 1 )) . 

However, at small separations, the situation becomes significantly more involved, since displace- 
ments of and interactions between topological defects lead to non-linear and many-body terms in 
the effective interaction potential. These short-range forces can be used to "synthesize colloidal 
molecules" [[SS] through a direct assembly using optical tweezers. Such aggregates are relatively 
stable objects made of dipolar and quadrupolar "colloidal atoms" which are bonded by topologi- 
cal defects. To obtain the free energy associated with a "chemical bond" , one needs to refine the 
description of the nematic mesophase and to take into account the behavior of topological defects 
at small separations. 

Several approaches which go beyond the linearized theory and a simple superposition approxi- 
mation have been proposed. Most popular are (i) numerical minimization of the Frank-Oseen elastic 
free energy functional in two and three spatial dimensions [ I34[ I71[ I72j (ii) minimization of the 
Landau-de Gennes free energy functional based on the five component tensorial order parameter [ 
E21 E31 E31 E31 EH EH ESI E3 EH] , (iii) molecular dynamics [[2SJ and Monte Carlo simulations [[32], 
or using classical density functional theory [[79] [80] . The results of these efforts can be summarized 
as follows. At small separations between particles the topological defects rearrange their spatial 
position in order to minimize the free energy. This structural rearrangement changes the character 
of the effective interactions as compared to the asymptotic multipolar expressions in Eqs. (|11[) . 
These changes are illustrated in Fig. [2] for the case of two long cylindrical and in Fig. [3] for two 
spherical colloidal particles. The mismatch in the boundary conditions at the surfaces of the parti- 
cles and the bounding box produces compensating topological defects: accompanying disclination 
lines in the case of cylinders, and equatorial rings in the case of spheres. The spatial arrangement 
of these defects depends on the distance between the particles. In both cases, the far-field director 
distribution has quadrupolar symmetry, which results in d~ A and d~ 5 decay of the interaction 
potential for cylinders and spheres, respectively. However, at some separation, which depends on 
the angle 9 and the particle type, repulsion between the particles changes to attraction. Switching 
between repulsion and attraction is more pronounced for cylinders than spheres, where it appears 
only for 9 = tt/2. The insets of Figs. [21 [3] show that such behavior is due to the rearrangement 
of topological defects. For spherical particles, for example, two Saturn rings (equatorial defects) 
merge into one at small separations. 
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dIR 

Figure 2. (Color online) Effective interaction potential (or potential of mean force) as a function 
of the separation d between two long cylindrical rods aligned along the z axis. Several angles 8 
between the far field director and the center-to-center vector are shown. The lower inset shows 
the z cross-section. The free energy is calculated by numerically minimizing the Landau-de 
Gennes free energy functional, Eqs. Q, ©. The far field director is along the x axis. The 
particle surfaces provide rigid homeotropic anchoring which gives rise to two disclination lines 
of topological charge of —1/2 at each particle. The alignment of the director and the variation 
of the scalar order parameter are shown in the top inset. The typical values for a nematic 
compound 5CB [ [70] are used in Eqs. @, ©: a = 0.06 x 10 6 J/m 3 , b = 0.816 x 10 6 J/m 3 , 
c = 0.45 x 10 6 J/m 3 , Li = 6 x KT 12 J/m, L 2 = 12 x HT 12 J/m, and L 3 = 0. K 2 = 9^ ulk L 1 /2. 

The nematic coherence length £ = (48Lic/6 2 ) ~ 10 nm, which corresponds to 5CB at a 
nematic-isotropic coexistence. The radius of colloidal particles is R/^ = 10. The system size is 
L x x L y = 40.R x 40R. 
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Figure 3. (Color online) Effective interaction potential between two spherical colloidal particles 
as a function of separation d. The geometry is sketched in the lower inset. As before, the Landau 
de-Gennes free energy functional, Eqs. Q, ((6|, is minimized numerically using the finite elements 
method and adaptive meshes. Particles provide rigid homeotropic boundary conditions for the 
director. The surface scalar order parameter is fixed to the bulk value of the order parameter, 
^surface = Sbuik- Material parameters used in Eqs. ((4J, © are similar to those used for cylindrical 
particles, see caption to Fig. [2] except that a = 0.01 x 10 6 J/m 3 . The upper inset depicts the 
nematic configurations for 9 — n/2 and several values of d. The blue rings around the particles 
show the iso-surfaces of the scalar order parameter corresponding to S = Sbuik/2. The radius of 
colloidal particles is R/£ — 5. The system size is L x x L y x L z = 30R x 30-R x 30i?. 
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5. Confinement effects and thin films 

One more example of a colloidal emulsion/ suspension with complex anisotropic interactions 
between inclusions is a system with a free standing smectic film serving as the host liquid. Here, 
the effective interactions can be controlled not only by anchoring conditions at surfaces of droplets 
or particles, but also by the symmetry of the mesophase, e. g. by using different smectic phases: 
SmA, SmC, or SmC*. Additionally, temperature, polarization, film thickness and droplet size can 
be used to tune these interactions [ 181] 182] 183] . In the following section, we briefly summarize the 
results on the self-assembly of inclusions in free standing smectic films. A more comprehensive 
overview of this field can be found in Ref . [ [HI] ■ 

5.1. Inclusions in free standing smectic films 

Thin free-standing smectic-C (SmC) or chiral smectic-C* (SmC*) films are ideal hosts for 
studying the assembly of inclusions in two-dimensional anisotropic colloidal suspensions and dis- 
persions. Indeed, there are several benefits in using thin films compared to bulk systems. First, a 
free-standing film is bounded by two liquid-air interfaces. This enhances the smectic order which 
may persist even at temperatures up to 30° C above the corresponding bulk SmC/SmC* - isotropic 
(or nematic, cholesteric) transition temperature and guarantees that the smectic ordering will not 
be destroyed upon heating of the film. Second, smectic ordering facilitates the interpretation of 
the optical film texture when observed between crossed polarizers. Finally, the free-standing film 
geometry allows one to use reflection as well as transmission optical modes for visual inspection of 
the films. 

Depending on a liquid crystal compound used, different types of inclusions can be nucleated in 
a film upon heating: (i) isotropic droplets close to the SmC-isotropic transition [ [55] , (ii) nematic 
(N) droplets nucleated in SmC films close to the bulk SmC-N transition [ [8TJ [86] |57| [88] , (iii) 
cholesteric (N*) droplets nucleated in a SmC* film close to the corresponding bulk phase transition [ 
I8T] [82] [83] 189] 150] . and (iv) smectic islands, or circular localized regions containing additional 
smectic layers [ 183 ] 188 ] 191] . 

Dispersions are normally prepared by using well-controlled cycles of heating/cooling of free 
standing smectic films [ 90 . In order to initiate the nucleation of lower order inclusions, the smectic 
film is gradually heated across the corresponding bulk transition temperature until the temperature 
of nucleation is reached. As soon as the nucleation process starts, the temperature is kept fixed. 
The process of nucleation depends on the film thickness. For thin films (Ni ayers < 24), initially 
nucleated droplets continue to grow by coalescing, which stops when the diameter of inclusions 
reaches the size of 4 — 5 /zm. For thicker films, nucleated droplets rarely coalesce, but continue to 
grow until the equilibrium sets in within the film, and droplets have grown to the diameter of 
about 18 — 20 /im [ 90 . Hence, the concentration and the size of inclusions can be controlled by the 
temperature difference, speed of a heating cycle, and by the initial film thickness. Upon cooling 
down, the inclusions can be transformed into smectic islands i. e. localized regions containing 
additional smectic layers [[55], 

Effective interactions between inclusions in a film depend strongly on the type of a liquid 
crystal. In a SmA phase, the n-director, which determines the average orientation of the long 
axis of molecules, is normal to the smectic layers. Hence, additional elastic deformations lead to 
long-range isotropic interactions between inclusions. In SmC phases, the n-director is tilted with 
respect to the layer normal by some angle a, as shown in Fig. The projection of the n-director 
onto the plane of smectic layers defines a two-dimensional vector, the c-director. Due to this tilt, 
SmC films provide an anisotropic environment for nucleated inclusions. The situation is even more 
complicated in SmC* films [[92]. This phase has the Ci symmetry with a two-fold rotation axis 
normal to the n-director and parallel to the smectic layers. Therefore, if a polarization has a 
component along the Ci axis, it can not be canceled by the two-fold rotations around the axis. 
This leads to a net polarization in the chiral SmC* phase in the p-direction as shown in Fig. 2}}. The 
coupling of this polarization to the c-director gives rise to a temperature-induced rearrangement of 
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Figure 4. (a) SmC and (b) SmC* phases. In the SmC phase the long molecular axis (or n- 
director) is tilted with respect to the layer normal at an angle a. SmC phase has monoclinic 
symmetry with the following group elements: (i) a plane of symmetry, the zy plane, (ii) a two-fold 
rotation axis perpendicular to this plane (a;-axis), and (iii) an inversion point O. A projection of 
the n-director onto the plane of smectic layers defines the c-director, which is a two dimensional 
vector. The SmC* phase is a chiral analogue of the SmC phase and is characterized by a helical 
rotation of the c-director. Due to the presence of the chiral molecules the c-director twists as one 
progresses through the layers. The pitch p is usually much larger then the inter-layer spacing. 
Chirality eliminates the symmetry plane as well as the inversion center and the remaining two- 
fold symmetry axis leads to a spontaneous polarization in its direction i. e. SmC* is a ferroelectric 
material [|92]. The vector p = e z x c defines the polarization direction. 



topological defect^] and can even change anchoring conditions at the boundaries of N* inclusions 
in free-standing ferroelectric SmC* films [ [521 183) . 

Smectic-inclusion interfaces impose specific anchoring conditions on the c-director. Isotropic or 
nematic inclusions in SmC films impose planar anchoring on the c-director [ [551 |BH1 [571 [Ml (55J ■ A 
similar situation is observed for smectic islands [IHH1ISI]- Cholesteric droplets nucleated in SmC* 
films impose homeotropic anchoring on the c-director [ 89, 90 , which can be rationalized by taking 
into account the spontaneous polarization of SmC*, which gives rise to a polarization-dependent 
term in the anchoring energy [ I94j . This term favors a planar orientation of the polarization or a 
homeotropic orientation of the c-director at the inclusion boundary, in order to achieve zero surface 
electric charge and to minimize the electrostatic energy of the inclusion. 

In some cases, anchoring conditions at the inclusion boundary can depend on temperature 
and thickness of the SmC* film or even a surface transition can occur [[52]. For example, during 
the nucleation process, just slightly above the bulk SmC*-N* transition, the c-director has planar 
anchoring. Upon further heating, a series of transitions from planar to homeotropic anchoring and 
back is observed until the film thickness reaches a critical value of N ss 9, below which only planar 
anchoring is possible. The transition from planar to homeotropic anchoring is initiated by splitting 
of the —1/2 surface defects into pairs of the —1/4 surface defects, which repel each other and slide 
along the inclusion boundary. Eventually, these intermediate defects of strength of —1/4 recombine 
either with the original intact —1/2 defect or with each other. 

A change in polarity of thin smectic films can also lead to a change of the type of defects 



1 Note that the accompanying topological defects of the c-director field are different from those of the n-director. 
Indeed, the c-director, contrary to the n-director, is a two-dimensional vector. Hence, its ground state manifold 
is a unit circle S , and the fundamental group is Z (integers). Therefore, in free standing smectic films, the type 
of a topological defect generated by inclusions with homeotropic or planar boundary conditions is either a bulk or 
surface point defect of strength of —1 or two surface defects of strength of —1/2. 
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(a) (b) (c) (d) 

Figure 5. Sketch of the c-director streamlines around a circular inclusion in a two-dimensional 
SmC* film. Planar anchoring is assumed at the boundary of the inclusion. The dipole- 
quadrupole-reversed dipole configurational transition is illustrated: (a) a single surface defect 
of strength of —1 forming a dipolar configuration of the c-director, (b) two surface defects 
of strength of —1/2 forming a mixed configuration with the decreased dipolar moment, (c) a 
quadrupolar configuration, and (d) a reverse mixed configuration made of two surface defects 
of strength of —1/2. This transformation can be induced by temperature variation and is re- 
versible upon heating and cooling of the film, assuming the temperature does not cross the 
thinning transition as the number of smectic layers decreases by one [ 1821183] . 

around the inclusions. For example, in a nonpolar membrane of SmC, the configuration of the c- 
director field is quadrupolar with two —1/2 surface defects, while it is dipolar with one —1 surface 
defect in the case of polar membrane of SmC* [HI]. A similar transformation between the dipolar 
and quadrupolar configurations can be performed reversibly by changing the temperature of the 
sample [ [HH H3] ■ The transformation of the c-director proceeds through splitting of the surface 
topological defect of strength of —1 into two surface defects of strength of —1/2 (see Fig. [5]). As 
the temperature changes, the —1/2 defects move along the boundary and may recombine on the 
opposite side of the inclusion resulting in the reorientation of the c-director at the boundary by 
an angle equal tt. 

Anchoring conditions determine the type and the topological charge of defects accompanying 
inclusions and therefore define their self-assembly behavior. For example, cholesteric droplets in 
SmC* films can provide homeotropic anchoring for the director. In this case, a formation of chains 
is observed which is due to the effective dipole-dipole interactions between the inclusions [ 1551150] . 
On the other hand, planar anchoring is observed for nematic inclusions in SmC films [ 86, 87]. In 
this case, the c-director field around each inclusion has a quadrupolar symmetry, which leads to the 
formation of kinked chains at low concentrations of inclusions and to two-dimensional structures 
with a square lattice at higher concentrations. In some situations, inclusions in SmC* films impose 
planar anchoring. In this case, for dipolar inclusions, chains orient perpendicular to the far field 
c-director. For a quadrupolar configuration, kinked chains and rhomboidal clusters of four particles 
are formed [HUES]. 

In spite of a significant amount of experimental work, theoretical understanding of these systems 
is rather limited. The main reason is the complexity of their description, which often requires 
techniques suitable to treat several length scales, from smectic layering to long-range distortions 
of the director field around inclusions. In addition, inclusions have soft boundaries, which pose 
additional problems when specifying boundary conditions. 

The effective elastic force between particles in free-standing SmC films was shown to have an 
asymptotic power- law of dr 1 decay, where d is the separation between the particles, and to be 
inversely proportional to the smectic-air surface tension 7 [|S5]. For small separations, a super- 
position approximation employing the analytical single-particle solution was used to calculate the 
effective interaction between inclusions with planar anchoring [|96]. For weak anchoring conditions, 
the global minimum of the free energy corresponds to a configuration where inclusions are in a 
contact and the center to center vector forms an angle of ~ 36° with the far field c-director. For 
anchoring strengths above a threshold value, there is an optimal separation between inclusions 
which increases with the anchoring strength and can reach a value of three inclusion diameters. 
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Figure 6. (Color online) Sketch of a quasi-two-dimensional nematic colloid sandwiched between 
two glass plates in a wedge- like geometry. The cell walls impose planar and the colloidal particles 
homeotropic anchoring on the director. In the thinner part of the cell colloidal particles generate 
quadrupolar configurations in the surrounding nematic field, while in the thicker part of the cell 
dipolar configurations are more stable. Hence, the dipole-quadrupole configurational transition 
is governed by the confinement. For ~ 2.32pm silica particles dispersed in 5CB, the threshold 
cell thickness is h c « 3.5/xm [ 1102] . 



Additionally, a pair of inclusions becomes more aligned with the far field c-director as the anchoring 
increases. 

The case of chiral SmC* can be modeled by introducing an additional elastic term, oc V x c, to 
the two-dimensional elastic free energy density of SmC [ 97, 98 . This term allows for a spontaneous 
bend deformation of the c-director and favors the dipolar defect configuration with two —1/2 
surface defects placed at one side of the inclusion. Without this term, the quadrupolar configuration 
of the c-director with two antipodal surface defects of strength of — 1 /2 was found to be energetically 
more stable [ [Ml [S3 DUD] ■ 

The aforementioned work treats inclusions as if they were rigid bodies, which is clearly an 
oversimplification for liquid-like inclusions, since the anchoring as well as the surface tension is 
coupled to the director orientation at the inclusion surface and hence to the structure of the 
surrounding film. The latter varies with temperature, droplet size, external field, etc. The shape 
of the soft inclusion in the SmC film with planar anchoring was studied using the two-dimensional 
elastic free energy supplemented by two phenomenological terms describing the anchoring energy 
and the line tension respectively [ 1101] . It was found that the anchoring term favors an elongated 
shape of inclusions, while the line tension term favors the circular one. For any finite anchoring 
strength there is a critical value of the line tension parameter at which the elliptical shape of the 
inclusion becomes unstable. 

5.2. Quasi two-dimensional nematic colloids 

A different route to control the two-dimensional crystallization of colloids is to use a confinement- 
induced change of the charge of a topological defect accompanying a colloidal particle [ [HU [Ml ESI 
l6"9l 1102] . To achieve this, a dispersion of particles is confined in a wedge- like cell. The surfaces of 
the cell provide planar anchoring, while particles impose homeotropic anchoring of the director, 
as shown in Fig. [5] As it turns out, the Saturn ring defect (quadrupolar symmetry) is stable in 
the thinner part of the cell, while the satellite defect (dipolar symmetry) is favored in its thicker 
part. At a threshold thickness, h c /R ps 3, colloidal particles with both dipolar and quadrupolar 
director configurations can coexist with each other [l69 j 1102] . Hence, the dipolar to quadrupolar 
transformation of the symmetry of the director field around the colloidal particles can be achieved 
by changing the cell thickness, not the particle size. 

For thin enough cells, where only quadrupolar colloidal particles are present, particles self- 
assemble into long kinked chains and only sometimes small crystallites of nine and less particles 
appear [[68]. Bigger crystals can be assembled by using laser tweezers. The resulting crystal struc- 
ture is shown in Fig. [7^. The measured strength of the "bonds" in such a crystal is rather weak 
compared to the case of dipolar colloids. For example, when a particle approaches the head of a 
kinked chain, the measured binding energy is of the order of lO^k^T for 2.32/im particles. For a 
lateral approach, the binding energy is of the order of 0.3 x 10 3 fcBT. Due to the anisotropy and 
weakness of the binding energy, the assembled crystallites are rather fragile and "melt" even for 
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Figure 7. (Color online) Sketch of two-dimensional colloidal crystallites composed of (a) 
quadrupolar [ [55] , (b) dipolar [ [ST] , and (c) mixed [ colloidal particles. The crystallites 
are assembled with the help of laser tweezers. 



small variations of temperature. 

Contrary to quadrupolar crystallites, the ones assembled from colloids with dipolar defects were 
found to be much more stable [[HI]- The measured binding energy for aligned dipoles is of the order 
of 5 x l0 3 kftT for 2.32//m particles (in both directions, 9 = and 7r/2). The arrays of long dipolar 
chains, forming a 2D colloidal crystal, were directly assembled as shown in Fig. [7|d. Each chain 
has a "ferroelectric" order of topological dipoles, and therefore the whole 2D assembly possesses 
"antiferroelectric" ordering. The effective interaction potential between two dipolar particles was 
also directly measured as a function of their separation d [ 61 . However, the accuracy of the 
measurements did not allow the comparison to the theoretical predictions, which concluded that 
the dipole-dipole power law of <i~ 3 is exponentially screened for distances larger than the thickness 
of the cell, with an extended crossover region between the two regimes [ 1103) . 

In the coexistence region, where both dipolar and quadrupolar configurations are stable, the 
corresponding trapping energy is of the order of 5 x 10 3 fcBT for 4fim particles [[55] |B5]. Close to 
a dipolar particle with the fixed orientation of its topological dipole, three different equilibrium 
configurations of the second quadrupolar particle are possible. This allowed the assembly of a 
variety of two-dimensional colloidal crystals with different unit cells, one of these crystals is shown 
in Fig. [7fc. The measured dipole-quadrupole pair interaction when particles approach each other 
along the rubbing direction (9 = 0) with the hyperbolic hedgehog defect placed between them, 
agreed well with the asymptotic behavior of d~ 5 . No measurable effects of the cell thickness on the 
power-law exponent was found, however, the prefactor of the interaction energy was found to be 
thickness dependent. 

The self-assembly idea was further extended to the field of "colloidal chemistry" [[69], where 
colloidal particles play the role of atoms, which could, by sharing their topological defects, bind 
to "colloidal molecules". Several types of such "colloidal molecules" were "synthesized" with the 
help of optical tweezers. A number of exotic structures, such as one-dimensional ordered wire- 
like structures were also assembled. The particles in such a wire are bound by a common closed 
disclination line, with the binding energy of the order of lO^k^T f° r 4.7/rai particles, which is 
comparable to the binding energy along a dipolar chain [ I104j . In addition, an externally applied 
AC electric field was used to obtain close-packed hexagonal lattices [ 1105] , 



6. Colloids at interfaces 



Liquid-liquid interfaces provide an additional way of controlling colloidal assembly. In case 
of a two-dimensional organization of colloidal particles at fluid interfaces, interparticle forces are 
normally isotropic and lead to a hexagonal ordering of particles. Anisotropic interactions are also 
possible and can be observed for non-spherical particles trapped by fluid- fluid interfaces [ 1106) . 
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In anisotropic fluids, such as liquid crystals, a much richer variety of ordered patterns at inter- 
faces can be observed, for example the chain-like assemblies of particles [ 11071 [TU51 1109j . The type 
and degree of ordering can be controlled by adding a molecular surfactant. For example, reorgani- 
zation of micro-particles was observed by the addition of an anionic surfactant SDS to the liquid 
crystal 5CB. First, the spacing between the particles within chains increased and ultimately two- 
dimensional arrays with local hexagonal symmetry were formed [ 1109] . These experiments demon- 
strated that reversible, chemosensitive control of the interfacial organization of colloidal particles 
is possible. An alternative method is to dope the emulsion with a photosensitive amphiphilic dye. 
Upon exposing the emulsion to ultraviolet irradiation, the dye undergoes trans-cis isomerization, 
which changes the nematic-air surface tension, anchoring on the surface of droplets, and elastic 
constants of the host. In this way one gains a control of the lattice constant of the two-dimensional 
hexagonal crystal formed by the glycerol droplets at the nematic-air interface [ IllOj . By irradiating 
the emulsion with spatially-modulated light one may even obtain patterned two-dimensional col- 
loidal crystals. More exotic situations include coexistence of different two-dimensional assemblies 
of glycerol droplets at a nematic-air interface, which can be controlled by reorientation of elastic 
dipoles around each droplet [ llllj . 

Apart from controlling the assembly of particles at interfaces, it is also possible to use nematic- 
isotropic (NI) interfaces for ordering colloids into three-dimensional structures. For example, one 
can obtain cellular solids by dispersing colloidal particles in the isotropic phase of a liquid crystal 
and then quenching it below the NI transition temperature [ II 121 [TT5] . Colloidal particles, densely 
packed at thin interfaces between different nematic domains, form an open cellular structure with 
a characteristic domain size of 10 — 100/^m. The size of the domains can be controlled by changing 
the particle concentration. Due to the cellular morphology, such gel-like solids have a remarkably 
high elastic modulus, which varies linearly with the particle concentration [ 1114) . 

The structural organization of colloidal particles can also be controlled by changing the speed 
of a moving NI interface. When a colloidal particle is captured by a NI interface, it can be dragged 
by it. As a result, periodic, stripe-like structures can be obtained, with the period depending on 
particle mass, size, and interface velocity [ 11151 ITTBl HIT) . In combination with externally applied 
electric field, a wide variety of particle structures is observed, ranging from a fine-grained cellular 
structure to stripes of varying periods, to a coarse-grained "root" structure [ 1118] . A combination 
of a moving nematic-isotropic interface with a patterned electric field can also be used to move 
particles from one place to another and pack them in a certain place in an ordered fashion. The 
speed of the interface, the magnitude of the applied electric field, particle size, density and its 
dielectric properties control which particles can be moved and which are left behind. This provides 
the possibility, for example, to place "defects" at particular locations in photonic crystals [ 1119) . 

From the theoretical standpoint, one can think of two origins of anisotropic interactions between 
particles pinned at an interface. First, the director deformations extend into the nematic phase, i. e. 
in addition to interfacial energies, the bulk elasticity contributes to the total surface free energy of 
the system. Second, the effective surface tension depends on the orientation of the director in its 
vicinity and hence varies along the interface due to the presence of colloidal particles. 

Long-range elastic effects significantly complicate theoretical understanding of interparticle 
interactions at NI interfaces: expressions for the free energy or force-distance profiles can hardly 
be obtained analytically, especially for small separations, the case when also many-body effects 
are important. Asymptotic analysis predicts that, for the case of the homeotropic anchoring at 
the interface, the elastic interaction is quadrupolar, repulsive, and decays with the distance as 
(R/d)~ 5 [ 11201 \TH\ . Numerical calculations are also not straightforward because the interfacial 
width and the size of the defect core are normally much smaller that the size of a colloidal particle, 
hence there are several rather different length-scales involved. Therefore, significant computational 
efforts are required in order to address either static or dynamic problems. The remedy to this 
problem is either to solve the discretized Euler-Lagrange equations or to directly minimize the 
free energy functional using finite elements with adaptive meshes [ \122\ [T23| 1124) . In what follows, 
we illustrate this approach by studying interaction of particles captured by a nematic-isotropic 
interface. 
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Figure 8. Schematic drawing of the geometry of the system considered. Nematic film of thickness 
h is sandwiched between the isotropic phase and a solid substrate, imposing unidirectional strong 
planar anchoring on the director. Weak homeotropic anchoring is assumed at the NI interface. A 
small temperature gradient is applied in the z direction, such that the free NI interface (without 
the particles) is stabilized at the z — plane (nematic is the thermodynamically stable phase for 
z < 0). Colloidal particles of radius R which are distance d apart are trapped at the interface. 
We consider both strong homeotropic and planar degenerate anchoring at particle surfaces. In 
the case of a single colloidal particle, we vary the ^-coordinate of the center of the particle. For 
a pair of particles, the z-coordinates of their centers are fixed at z — 0, and the separation d is 
varied along the x direction. 

6.1. Example: Colloidal particles at a nematic-isotropic interface 

In this section, we illustrate how to determine the force on a colloidal particle from the side 
of a NI interface, as well as evaluate the pair force between two colloidal particles pinned at this 
interface. In order to avoid the divergence of the elastic free energy in the defect core, we use the 
expression for the free energy density in the Landau-de Gennes form, Eqs. Q, ©• This free energy 
has to be minimized either under constraints imposed by the inclusion (colloidal particle) or by 
adding an appropriate surface term to the total free energy. A similar two-dimensional problem 
has already been treated in Refs. [ 122, 1123] . In what follows, we consider the three-dimensional 
case, with a flat unperturbed NI interface and spherical colloidfl The geometry is sketched in 
Fig. [5] and is similar to the one studied experimentally in Refs. [ 11071 1111) . We assume hybrid 
boundary conditions: a strong parallel anchoring at the substrate and a weak homeotropic one at 
the NI interface. The latter is tuned by choosing the appropriate anisotropy of elastic constants, 
L2/L1 = —1/2 [ 1125] . We first calculate the force on a single particle as a function of its immersion 
in the nematic film. Two types of anchoring at the particle surfaces, homeotropic and planar, are 
considered, as well as different thicknesses of the nematic film. The planar degenerate anchoring 
at particle surfaces is modeled by the expression proposed in Ref. [ 126 . We then proceed by 
calculating the effective interaction potential between two particles pinned at the interface and 
discuss bridging forces and defect configurations around the colloidal particles. To summarize, we 
compare our results with theoretical predictions at large separations as well as to experiments. 

Single particle at a NI interface 

Let us first address the force on a single particle due to bending of a nematic-isotropic interface 
and elastic director deformations inside the nematic layer. The results of the minimization of 
the Landau-de Gennes free energy functional as given by Eqs. flU, © are summarized in Fig. |9] 
Figs.[9K,c show several cross-sections of the scalar order parameter S, which is the largest eigenvalue 
of the order tensor S a p. As one can see, the structure of the topological defects is different for two 

2 Note that the interaction between long rods considered previously in Refs. [ 28, 122 123 124 differs quanti- 
tatively from the interaction between spherical particles because another type of topological defects accompanies 
spherical colloids: either a point defect of strength -1 or a ring defect of strength -1/2 (contrary to -1/2 strength 
line defects for cylinders). 
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Figure 9. Contour plots of the nematic order parameter for (a) planar and (c) homeotropic 
anchoring of the nematic director on the particle surface. The particle position is fixed at a 
distance z/R = 3/5 above the NI interface. Plots (b) and (d) show the free energy as a function 
of particle immersion into the nematic film for planar and homeotropic anchoring, respectively. 
Several film thicknesses are shown. Insets illustrate the force on the particle evaluated by nu- 
merical differentiation of the corresponding free energy at z — as a function of the thickness 
of the nematic film. The twist elastic constant, K2 = 9S^ ulk Li/2, is evaluated at the bulk NI 
transition temperature. The following values for the material parameters are used in Eqs. @, 
©: b = 0.816 x 10 6 J/m 3 , c = 0.45 x 10 6 J/m 3 , L x = 6 x 10~ 12 J/m, L 2 = -3 x 10~ 12 J/m, 
Li = 0, and a — ajvi(l + z x 0.004), where clni = 6 2 /(24c) corresponds to the nematic-isotropic 
coexistence. In this case, R/£ = 5, where £ is the nematic coherence length defined in the caption 
to Fig. [2 The system size is L x x L y = 307? x 307?. 
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types of anchoring. In the case of homeotropic anchoring, a half Saturn ring defect is created in 
the nematic layer (the second half is merged into the isotropic phase), Fig. while for planar 
anchoring the defect structure is less symmetric and the defect core is split between two branches, 
each of them representing a disclination line which starts at the particle surface and ends in the 
isotropic phase. This defect cannot be easily identified with the standard structures, such as a 
Saturn ring, a satellite, or a bojoom-like topological defects. 

Different defect structures, however, do not lead to a quantitatively different dependence of 
the free energy on the particle immersion depth, which is shown in Fig. [5jD,d. In both cases, the 
particle is expelled from the nematic phase with a force that depends on the films thickness h. For 
small immersion depths, the free energy is roughly quadratic in the depth i. e. the force on the 
particle is proportional to the depth. The expelling force is also larger for homeotropic anchoring, 
since this type of anchoring introduces more deformations into the nematic phase. 

The asymptotic analysis for similar hybrid boundary conditions (both assumed to be rigid) 
predicted a power-law decay of the expelling force of at least (R/h) 8 [ 1121] . which basically reflects 
the fact that stronger elastic deformations of the director field are created in thinner nematic 
films. It is clear from the insets of Fig. [9] that this scaling does not hold in our case. The index 
of the exponent changes by practically six orders of magnitude. Note, however, that this might 
be an artifact of the assumption of the infinitely strong director anchoring at the particle surface, 
which results in stronger elastic deformations for thinner films. In addition, our calculations assume 
that the homeotropic anchoring at the NI interface is rather weak, contrary to infinitely strong 
homeotropic anchoring assumed in Ref . [ II 21] . 

Interaction of two particles pinned at a NI interface 

We now look at the effective interaction potential of two particles pinned at a NI interface. 
We adopt the same geometry as for a single particle, with both particles fixed at z = and a 
separation d. The corresponding geometry is sketched in Fig. El 

The first observation to make is that the structure of the director field, and hence the topological 
charge of the defects, is similar to the one observed for a single colloidal particle, at least at large 
separations. The free energy, however, has a very different dependence on the particle-particle 
separation. In both cases, it is practically insensitive to the film thickness h. For homeotropic 
anchoring, Fig. HUH , particles always repel each other. For finite h the director field around a 
single particle is given asymptotically by the dipole term [ I121j . which leads to a dipolar, (R/d) 3 , 
decay of the effective interaction at large distances. In our case, the fit of a power law function to 
the free energy dependencies predicts (qualitatively) a (R/d) 4 decay of the interaction potentials, 
indicating that the situation is more complicated than the theoretical prediction. 

A completely different situation is observed for planar anchoring, which is shown in Fig. llOb . In 
this case, particles weakly attract each other with a force which is almost independent of separation. 
However, below a threshold distance d c w 3i? which depends on h, this weak attraction increases 
by roughly a factor of ten. In fact, a similar attraction has already been reported earlier and is 
due to the formation of a nematic bridge or a bridging defect between the particles [ 11271 1123] . It 
effectively leads to an additional force due to the surface tension of the bridge. Capillary bridging 
leads to a hysteresis in the force-distance profile at small distances [ 1123] and should not affect the 
long-range interaction of colloidal particles. 

In conclusion, within the Landau-de Gennes formalism we have studied the effective pair inter- 
action between two spherical particles trapped at a NI interface. The interparticle force depends on 
the type of anchoring as well as the formation, annihilation, and interaction of defects. Therefore, 
one expect that many-body contributions will be important to properly describe collective colloidal 
ordering at interfaces [ 11241 1128] . 

7. Summary and outlook 

We have discussed theoretical, computational, and experimental techniques used to study var- 
ious inclusions in liquid crystals. In particular, special attention has been paid to two-dimensional 
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Figure 10. Contour plots of the nematic order parameter for (a) planar and (c) homeotropic 
anchoring of the nematic director on the particle surfaces. The particles are fixed at z = and 
separated by distance d/R = 2.7. The corresponding geometry is depicted in Fig. [8] Plots (b) and 
(d) show free energy as a function of distance d between the particles for planar and homeotropic 
anchoring, respectively. Several film thicknesses are shown. The twist elastic constant, K2 = 
9S^, ulk Li/2, is evaluated at the bulk NI transition temperature. The values of the material 
parameters used in Eqs. Q, (0 are given in the caption of Fig. [5] 
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systems, such as free-standing smectic films, thin planar nematic cells, and nematic-isotropic in- 
terfaces. 

From a computational point of view, we find that the Landau-de Gennes free energy functional 
based on a five component tensorial order parameter is the most suitable description of colloidal 
systems with anisotropic hosts. By construction, it removes singularities of the director field in 
the cores of disclination lines or point defects by allowing biaxial perturbations [ I129j and local 
melting of the mesophase. It predicts the internal structure of topological defects, is appropriate for 
describing nematic-isotropic interfaces and suitable for calculating short-range effective interactions 
between inclusions. The model has been used extensively to resolve defect structures accompanying 
a pair of colloidal particles at small separations, to calculate effective interaction potentials and 
associated binding energies between inclusions, and to study self-assembly of inclusions at liquid- 
liquid interfaces. 

In conclusion, the following issues still remain unresolved: (i) understanding of the temperature 
dependence and interplay of polarization and elasticity in SmC* films, leading to rearrangements 
of surface defects and changes in the type of anchoring at inclusion boundaries; (ii) construction of 
effective Hamiltonians for interfaces, which can account for bulk degrees of freedom in an effective 
way to avoid costly calculations using adaptive meshes; (iii) development of adaptive multiscale 
approaches, where classical density functional or molecular dynamics simulations are used on small 
scales, for example at the defect core, while a continuum model is used to describe slow variations of 
an order parameter on larger scales; (iv) development of bottom-up descriptions, which start from 
a chemical structure and atomistic force-field parametrization and continue with the construction 
of a coarse-grained particle-based representation [ 1130) , which is then used to parameterize an 
appropriate free energy functional. 
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